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SECTION  1 


INTRODUCTION 


A  space  based  radar  (SBR)  that  must  search  for  and  track  targets 
moving  close  to  the  earth's  surface  must  be  designed  to  discriminate  be¬ 
tween  these  targets  and  radar  clutter  caused  by  the  undesired  signal 
returned  from  the  earth.  For  targets  that  are  very  close  to  the  earth, 
clutter  and  target  returns  can  appear  at  the  same  range  in  the  radar  re¬ 
ceiver.  If  this  condition  occurs,  an  MTI  or  a  pulse  doppler  radar  may  be 
used  to  discriminate  on  the  basis  of  doppler  measurements.  However,  since 
the  clutter  return  usually  exceeds  the  target  return  by  many  orders  of 
magnitude,  care  must  be  taken  during  the  coherent  processing  to  avoid  ef¬ 
fects  such  as  aliasing  and  spectral  leakage.  Thus  an  accurate  simulation 
of  a  space  based  radar  must  include  an  accurate  model  of  land  and  sea 
clutter  in  order  to  evaluate  the  performance  of  different  coherent  proces¬ 
sing  techniques. 

The  purpose  of  this  report  is  to  describe  a  clutter  model  inten¬ 
ded  for  implementation  in  a  pulse  doppler  tracking  radar  simulation.  A 
standard  approach  to  clutter  modeling  (e.g.  Reference  1)  involves  division 
of  the  ground  area  illuminated  by  the  radar  antenna  into  a  number  of 
cells.  A  discrete  scatterer  of  known  cross  section,  amplitude,  doppler 
frequency  and  location  is  assigned  to  each  cell  to  represent  all  returns 

from  a  differential  area  element.  These  discrete  scatterers  or  clutter 

targets  then  enter  a  receiver  simulation  in  the  same  way  as  returns  from 
any  other  target  (e.g.  an  aircraft).  Advantages  of  this  approach  are  that 

it  is  st rai ght f orward  to  implement  and  that  it.  can  be  made  increasingly 

accurate  by  using  a  greater  number  of  cells.  Unfortunately,  because  of 
the  large  area  illuminated  by  an  SBR,  the  large  doppler  extent,  and  the 


resolution  required,  an  unacceptably  high  number  of  clutter  targets  would 
be  required. 

Here  a  second  approach  has  been  adopted.  This  technique  con¬ 
sists  of  modeling  the  clutter  as  a  random  process  with  mean  power,  mean 
doppler,  and  doppler  spread  determined  from  experimental  measurements  and 
from  physical  and  geometrical  quantities  pertinent  to  the  radar  location, 
look  angle,  beam  pattern,  etc.  The  above  three  quantities  are  used  to  de¬ 
fine  the  clutter  spectrum  observed  at  a  given  radar  range  gate,  and  a  sta¬ 
tistical  signal  generation  technique  is  employed  to  generate  a  realization 
or  sample  function  of  the  instantaneous  clutter  signal  observed.  This 
method  permits  simple  interpretation  of  results  and  reduces  simulation 
running  time  since  only  one  target  requires  processing.  The  accuracy  can 
be  increased  with  the  level  of  detail  of  the  analytic  clutter  model. 


SECTION  2 


MOOElING  APPROACH 


The  goal  of  this  work  is  to  provide  a  clutter  model  appropriate 
to  a  high  resolution  pulse  doppler  SBR.  It  is  desired  that  the  model 

provide  realizations  or  sample  functions  of  the  clutter  video;  that  is, 
the  model  should  provide  the  clutter  voltage  and  phase  as  a  function  of 
range  gate  for  all  the  pulses  of  a  coherent  pulse  train.  It  is  convenient 
to  separate  the  clutter  signal  returned  through  the  main  beam  of  the 
antenna  from  that  returned  through  the  antenna  sidelobes.  The  footprint 
on  the  earth's  surface  defined  by  the  3  dB  contour  of  the  two-way  antenna 
pattern  defines  the  area  that  contributes  to  the  main  beam  signal;  the 
remaining  area  contributes  to  the  sidelobe  signal.  Of  course  the 
illuminated  areas  are  further  limited  by  range  gating,  and  ambiguous  range 
returns  must  be  accounted  for.  These  effects  are  discussed  in  Section  4. 

Because  of  the  fortunate  statistics  of  the  clutter  signal 

observed  by  the  radar,  it  is  possible  to  utilize  a  simple  Fourier 

transform  technique  to  generate  realizations  of  the  clutter  video  signal 
observed  at  a  range  gate.  Given  an  analytical  or  numerical  description  of 
the  mean  clutter  spectrum,  in  this  technique  a  realization  or  sample 
function  of  the  clutter  voltage  is  obtained  by  taking  the  discrete  Fourier 
transform  of  a  specified  random  sequence.  The  necessary  mathematics  and 
justification  of  this  method  are  fully  described  in  the  Appendix.  An 
essential  assumption  is  that  the  quadrature  components  of  the  clutter 

voltage  are  uncorrelated  Gaussian  variates. 

In  order  to  simplify  the  sea  clutter  model  developed  herein, 
it  is  assumed  that  the  clutter  power  spectrun  is  well  represented  as  a 


Gaussian  in  frequency  whose  power  and  spectral  spread  may  be  computed  on 
the  basis  of  the  S8R -earth  geometry  with  consideration  of  the  illumination 
and  receiving  pattern  of  the  SBR  antenna. 

The  clutter  return  has  been  found  to  be  uncorrelated  from  one 
burst  waveform  to  another  when  the  RF  carrier  frequency  is  changed  by  at 
least  the  reciprical  of  the  pulse  width2.  Since  the  normal  mode  of  opera¬ 
tion  of  an  SBR  is  to  change  radar  frequency  between  pulse  trains,  this  de- 
correlation  property  is  retained  here  in  the  generation  of  the  clutter  for 
successive  pulse  trains. 

In  practice,  the  mean  clutter  spectrum  is  obtained  in  the  simu¬ 
lation  at  the  start  of  each  dwell  or  coherent  pulse  train.  It  is  assumed 
that  the  geometry  does  not  change  during  the  coherent  integration  time 
(dwell  duration)  so  that  the  clutter  spectrum  is  constant  during  this 
period.  For  each  dwell  a  sequence  of  numbers  is  drawn  from  a  random  num¬ 
ber  generator  and  a  realization  of  the  clutter  signal  is  generated  as 
described  in  the  Appendix.  Succeeding  dwells  are  assumed  to  utilize  a 
different  radar  frequency  and  thus  require  a  new  sequence  of  random  num¬ 
bers  to  generate  independent  clutter  signals.  Thus,  for  each  coherent 
pulse  train,  the  clutter  signal  at  a  particular  range  bin  is  determined 
using  a  statistical  signal  generation  technique. 

The  above  discussion  applies  to  clutter  observed  through  the  an¬ 
tenna  main  beam  as  well  as  to  clutter  that  enters  the  receiver  through  the 
antenna  sidelobes.  In  the  case  of  sidelobe  clutter,  it  will  he  shown  that 
a  flat  power  spectrum  is  appropriate  and  only  the  clutter  power  is  re¬ 
quired.  In  addition,  simpler  signal  generation  techniques  are  feasible 
that  do  not  require  Fourier  tranforms. 

Our  discussion  thus  far  has  focused  on  modeling  the  signal  in 
a  single  range  bin.  When  range  tracking  is  performed,  the  signals  in 
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several  adjacent  bins  must  be  modeled  and  it  is  important  to  include  the 
effect  of  bin-to-bin  correlation.  For  sea  clutter,  experimental  data 
reported  over  a  wide  ranqe  of  pulse  widths  and  carrier  frequencies  demon¬ 
strates  that  the  clutter  returns  decorrelate  over  a  ranqe  displacement  of 
about  one  pulse  length,  when  the  range  is  measured  in  units  of  delay  . 
This  is  understood  by  realizing  that  since  many  independent  scatterers 
contribute  to  sea  clutter  returns,  the  correlation  for  two  returns  separ¬ 
ated  in  ranqe  is  approximately  equal  to  the  fraction  of  illuminated  ground 
area  which  is  shared  by  the  returns.  For  rectanqular  pulses  with  low 
sidelobes,  the  correlation  is  triangular  and  qoes  to  zero  at  a  ranqe 
displacement  equal  to  a  pulse  length.  Thus  the  sea  clutter  returns  will 
be  uncorrelated  for  ranqe  bins  separated  by  more  than  a  (compressed)  pulse 
length. 


In  qeneral,  the  effects  of  ambient  refraction  are  unimportant 
for  the  clutter  model  described  here.  Even  for  low  altitude  orbits  the 
path  length  over  which  refraction  occurs  comprises  a  very  small  fraction 
of  the  total  satellite  to  earth  path.  Thus,  the  fractional  change  in 
radar  range  is  very  small  and  will  not  s iqn i f icant ly  impact  calculation  of 
the  illuminated  areas  or  the  radar  ranqe  equation.  The  larqest  effect  is 
an  increase  of  clutter  reflectivities  that  occurs  with  increasing  grazinq 
anqle.  This  effect  is  important  only  for  qrazing  angles  less  than  a  few 
degrees;  fortunately  this  ranqe  is  outside  that  currently  under  considera¬ 
tion  for  SRR  operation.  [f  S3R  operation  at  very  small  grazinq  anqles  is 
to  be  included  here,  the  calculation  of  qrazing  angles  must  be  modified  to 
include  deviations  caused  by  refraction. 

Sianal  doopler  effects  are  dependent  on  the  earth's  surface 
rotational  velocity  and  this  has  been  included  in  our  f ormul at ’on.  First, 
for  the  sake  of  simplicity,  we  treat  the  doppler  effects  under  the  assumD- 
tion  of  a  non-rotating  earth.  Later,  the  approach  is  extended  to  the  case 
of  a  rotating  earth.  In  this  reqarq  the  radar  look  anqles  are  measured 
relative  to  the  satellite's  velocity  vector  which  lies  in  the  horizontal 
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plane  as  shown  m  Figure  1.  In  conformance  with  the  later  discussion 
which  includes  the  effects  of  rotation  of  the  earth  and  of  vertical 
satellite  motion,  the  elevation  angle  (EL)  is  negative  for  a  downward 
looking  1 ine-of-s ight  (LOS). 

A  pulse  doppler  radar  is  able  to  mitigate  the  effects  of  clutter 
by  separating  target  and  c’utter  in  the  dopple’"  freguency  domain  through 
the  use  of  a  large  number  of  coherent  pulses  repeated  in  tine  at  a  high 
pulse  repetition  frequency  (PRF).  An  incteosc  in  the  PRf  gives  an 
increase  in  the  total  doppler  bandwidth  surveyed.  /n  increase  in  the 
number  of  pulses  increases  the  power  as  well  as  the  doppler  resolution. 
Pulse  amplitude  weighting  on  receive  is  utilized  to  suppress  doppler  side- 
lobes  associated  with  the  coherent  processing.  The  effect  of  a  particular 
clutter  target  depends  upon  the  radar  velocity,  wavelength  and  the  angle 
of  the  clutter  from  the  antenna  boresight.  Since  returns  from  clutter  are 
weighted  by  the  antenna  pattern,  improvement  in  clutter  rejection  can  be 
obtained  by  narrowing  the  antenna  beam.  This  technique  not  only  realizes 
less  overall  clutter  power  but  also  decreases  the  spect-al  spread  of  the 
clutter  contribution  by  selectively  renoving  returns  from  high  off-axis 
angles. 

A  quadratic  phase  shift  applied  on  a  pulse-to-pulse  basis  can  be 
implemented  in  order  to  mitiqate  the  spread  in  the  main  beam  clutter 
signal  due  to  ambiquous  range  returns  (Reference  2).  Quadratic  phase 
shift  processing  (QPSP)  ’s  discussed  in  detail  in  Section  4.  Improvement 
in  clutter  rejection  can  be  large  and  depends  upon  the  SSR  look  angle  and 
the  relative  importance  of  platform  motion  on  the  clutter  spectral  width. 

Another  technique  for  reducinq  platform  motion  caused  doppler 
spread  which  is  not  investigated  here  is  the  use  of  displaced  phase  center 
antennas  (Skolnik,  Reference  4;  Knepp  and  Dana,  Reference  5). 
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Figure  1.  Clutter  calculation  geometry. 


SEA  CLUTTER  REFLECTIVITY  MODEL 


Sea  clutter  is  the  radar  echo  received  as  a  result  of  reflec¬ 
tions  from  the  sea's  surface.  The  effective  radar  cross  section  (RCS)  of 
the  sea  echo  may  be  predicted  by  taking  the  product  of  the  sea  reflectivi¬ 
ty  and  the  area  that  contributes  to  the  echo  strength.  Sea  clutter  re¬ 
flectivity  is  usually  reported  in  terms  of  the  parameter  o0,  the  mean 
radar  cross  section  per  unit  area  of  illuminated  surface.  Thus,  a  value 
of  0q  =  -30  dB  corresponds  to  an  average  reflectivity  of  lu-3.  o0  is  a 
function  of  grazing  angle  (defined  as  the  angle  between  the  tangent  to  the 
surface  and  tne  radar  LOS),  polarization,  frequency,  and  sea  state. 

Because  of  the  wide  range  of  SBR  look  angles  of  interest,  a  mo¬ 
del  of  Oo  over  the  range  of  grazing  angles  from  a  few  degrees  up  to  normal 
incidence  is  required.  Data  at  high  grazing  angle  is  necessary  to  esti¬ 
mate  the  sidelobe  altitude  return,  although  the  SBR  main  beam  would  be 
overwhelmed  by  clutter  at  grazing  angles  greater  than  50  or  60  degrees. 
The  model  we  have  determined  for  og  is  shown  in  Figure  2  as  a  function  of 
grazing  angle  at  L-band  (1.5  GHz).  Two  curves  are  shewn,  the  high  wind 
(HW)  and  low  wind  (LW)  models,  with  the  difference  in  the  two  representing 
the  effect  of  heavy  and  light  seas,  respectively.  The  equations  for  this 
model  are 

o0  =  7. 0-10-  V*01  +  0.841  exp  | -(  x/2-*)2/ .0247  }  ;  HW  (1) 

o0  =  3.6*  10- V* 72  ♦  2.810  exp  j-(  "/2-iT)7/.0149}  ;  LW  (2) 

where  I  is  the  grazing  angle  in  radians.  The  first  terms  in  the  above  are 

the  high  wind,  low  wind  models  of  Tomlinson6  evaluated  at  L-band.  We  have 


found  them  to  provide  a  reasonable  fit  to  reported  data3*7  averaged  over 
polarization  and  for  grazing  angles  under  about  seventy  degrees.  The 
additional  Gaussian  terms  have  been  added  to  provide  a  closer  fit  to  data 
at  high  grazing  angles. 


SECTION  4 


MAIN  BEAM  CLUTTER  MOOEL 


Implementation  of  the  method  outlined  above  requires  calculation 
of  the  mean  power  and  doppler  spectrum  of  the  clutter  signal.  Temporarily 
neglecting  the  effect  of  range  ambiguities,  the  Illuminated  area  depicted 
In  Figure  1  Is  defined  in  azimuth  by  the  beamwldth  and  In  elevation  by  the 
smaller  of  the  elevation  beamwldth  or  pulse  length.  These  angular 
variations  cause  a  doppler  spread  In  the  clutter  signal  through  coupling 
to  the  SBR  platform  motion  In  the  doppler  equation 


2vs 

f  *  — -  cos (EL)  cos(AZ) 


(3) 


2v 


Af ^  cos(EL)  sln(AZ)  AAZ 


(4) 


2vs 

AfEL  *  “x"  sin  (EL)  cos(AZ)  aEL  . 


(5) 


Equation  3  gives  the  mean  doppler  frequency  as  a  function  of  the  antenna 
geometry.  Equations  4  and  5  give  the  3dB  spectral  widths  In  azimuth  and 
elevation,  respectively.  In  the  above  equations 


AAZ  =  3  dB  two-way  azimuth  beamwidth  (AZ3) 
v$  *  satellite  velocity 


and 


&EL  =  minimum  (EL5,ct  tan(1/)/2K) 


(6) 


where 


EL 3  =  3  dB  two-way  elevation  beanwidth 
c  =  speed  of  light 
tc  =  compressed  pulse  length 

R  *  slant  range 


Because  of  the  high  satellite  velocities,  the  coupling  between 
beamwidth  and  platform  motion  is  usually  the  most  Important  effect  in 
determining  doppler  bandwidth.  Since  the  antenna  gain  may  be  approximated 
by  a  Gaussian  dependence,  the  azimuth  and  elevation  doppler  spreads  are 
also  Gaussian  and  their  variances  may  be  summed  to  yield  a  net  spectral 
spread  due  to  beamwidth  effects 


o 


B 


2 


(7) 


where  the  standard  deviations  are  related  to  the  3  dB  spectral  widths  by 


'AZ  ^  *42 

«AZ 

(8) 

’EL  =  *42 

a,el  • 

(9) 

Internal  or  wavp  motion  of  the  sea  also  broadens  the  doppler 
spectrum.  Many  measurements  of  sea  clutter  have  yielded  a  Gaussian  spec¬ 
trum  whose  spread  is  related  to  the  internal  velocity  spectral  width  of 
spa  clutter,  o^,  as 


o 


l 


(10) 


This  author  is  not  aware  of  a  good  data  base  to 
that  ar?  applicable  to  the  large  illuminated 
Existing  measurements  of  the  mean  doppler  shift 
ize  the  smaller  resolution  cells  of  airborne 


predict  values  of  oy 
areas  viewed  by  an  SBR . 
and  spectral  spread  util- 
and  ground  based  radars. 
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Figures  3  and  4  show  data  for  the  spread,  a  ,  and  mean  shifts  respectively 
as  a  function  of  sea  state.  Typical  high  wind  values  for  the  spread  and 
mean  velocities  are  seen  to  be  about  1  m/sec  and  3  m/sec,  respectively. 
Roughly  speaking,  for  an  SBR  the  mean  shifts  will  be  converted  to  a  spec¬ 
trum  spread  through  the  integration  of  a  large  number  of  smaller  sized 
cells.  Thus  we  can  estimate  an  Internal  velocity  spread  of  3  or  4  m/sec 
for  sea  clutter  viewed  by  an  SBR.  These  numbers  are  consistent  with  the 
value  of  5  m/sec  suggested  by  Tomlinson6.  The  accuracy  of  this  estimate 
will  be  improved  as  SBR  measurement  data  becomes  available. 

Combining  the  effects  of  platform  motion  and  internal  motion 
yields  a  net  doppler  variance  given  by 


We  now  consider  the  effects  of  ambiguous  range  returns.  Figure 
5a  depicts  the  illuminated  areas  associated  with  the  ambiguous  range 
returns  for  a  particular  range  bin.  Doppler  spectra  are  associated  with 
the  signal  power  returned  from  each  one  of  these  ambiguous  clutter  returns 
as  discussed  above.  The  signal  voltage  for  this  range  bin  is  the  sum  of 
the  signal  voltages  from  the  ambiguous  returns.  In  the  following  we  will 
refer  to  this  simply  as  the  sum  voltage  and  to  its  doppler  spectrum  as  the 
sum  spectrum.  Rather  than  generate  statistical  signal  samples  represent¬ 
ing  each  ambiguous  return  and  perform  this  sum  directly,  here  we  analyti¬ 
cally  compute  the  sum  spectrum  so  that  only  the  sum  voltage  requires 
statistical  sampling  on  a  pulse-to-pulse  basis. 


'M 


Returning  again  to  the  individual  spectra  of  the  ambiguous 
returns,  we  note  that,  because  of  the  dependence  of  elevation  angle  on 
range  and  also  the  dpper, dence  of  doppler  frequency  on  elevation  angle,  the 
mean  doppler  or  center  frequency  of  these  spectra  will  vary  over  the 


_t 
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Measured  spectral  spread  of  sea  clutter 


ambiguous  returns.  There  exists  a  signal  processing  technique,  QPSP,  that 
can  eliminate  this  doppler  shift  between  ambiguous  range  returns.  An 
explanation  of  how  QPSP  can  be  implemented  will  be  given  at  the  end  of 
this  section.  For  now  we  assume  that  QPSP  can  mitigate  the  doppler  shifts 
between  ambiguous  returns  as  stated,  and  continue  to  investigate  the 
properties  of  the  sum  signal. 

Without  QPSP  the  doppler  spectral  width  of  the  sum  siqnal  due  to 
ambiguous  returns  is  determined  by  the  antenna  elevation  beamwidth,  since 
the  center  frequencies  of  the  ambiguous  returns  vary  with  elevation 
angle.  Since  the  individual  spectra  are  independent,  the  net  power  is  the 
sum  of  the  individual  powers  of  all  the  ambiguous  returns. 

With  optimum  implementation  of  QPSP,  the  spectra  of  the  ambigu¬ 
ous  returns  have  identical  center  frequencies  as  shown  in  Figure  5b.  For 
this  case  the  net  power  is  still  the  sum  of  the  individual  powers,  but  the 
doppler  spread  of  the  sum  signal  is  the  same  as  that  computed  for  one 
ambiguous  return.  This  may  be  seen  in  the  following. 

The  autocorrelation  of  the  voltaqe  from  the  k'th  anbiguous 
return  may  he  expressed  as 


where 


\(t)  *  PkR(  c) 

=  '*k(t)  *kU+  0> 

=  moan  power  of  the  k't'n  return 
*k  =  voltage  of  toe  k’tn  return 


(1?) 

(13) 


This  formulation  is  an  explicit  statement  that  the  shape  of  the  individual 
spectra  are  identical  and  lienee  the  shape  of  the  autocorrelation  functions 
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are  also  Identical  since  the  antenna  azimuth  beamwidth  is  approximately 
constant  over  the  antenna  main  lobe.  The  P^'s  are  not  all  equal  since 
they  depend  on  the  antenna  elevation  gain  which  varies  with  ambiguous 
range.  The  autocorrelation  of  the  sun  siqnal  from  the  combined  ambiquous 
returns  is  then 

R£(t)  .  <(E  Xk(t))(E  x*,(t+T))>  (14) 

k  k 

=  <E  xk(t)  x*  (t+t)>  (15) 

k 

since  the  ambiguous  returns  are  independent.  Making  use  of  Equations  12 
and  13,  the  autocorrelation  of  t'ne  sum  voltage  with  QPSP  may  be  expressed 
as 

Rz(x)  =  r(t)Z  Pk  .  (16) 

Since  the  shape  of  the  autocorrelation  function  qiven  by  Equation  16  is 
the  same  as  that  of  the  individual  autocorrel ation  function,  the  shape  of 
the  sum  spectrum  will  be  the  same  as  that  of  the  individual  spectrum. 
Therefore  our  modeling  assumption  is  verified. 

We  now  list  the  algorithms  and  equations  needed  to  implement 
the  above  formulation. 

4.1  MAIN  BEAM  CLUTTER  POWER  IN  THE  M'TH  RANGE  BIN 

In  this  subsection  an  expression  is  given  for  the  clutter 
contribution  to  the  power  received  in  the  m'th  range  bin  from  all  clutter 

sources  within  the  footprint  of  the  main  beam.  Several  prerequisite 
calculations,  are  first  performed. 
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The  range  to  the  center  of  the  main  beam  footprint  is  given  by 
=  .  (r0+h)  sin(EL0)  -  /(r  Q+h)  2  sin  2(F.L  0)  -  (2r0h+h*)  (17) 

where 

r0  =  earth's  radius 
h  =  satellite  altitude 

EL0  =  elevation  angle  of  the  LOS  . 

The  grazing  angle  at  the  footprint  center  is 

V  =  sin- 1  jjl  (l+h/2r0)  -  -^j  <i8) 

jRc  2r o) 

and  the  chanqe  in  range  over  the  elevation  beamwidth  is 

AREL  =  ?Rc  si"(f-L  3/2)1  /  {cos(EL3)-  cos  (2*)}  (10) 


The  minimum  and  maximum  ranges  to  the  footprint  are 


R  . 

min 

R 

max 


AR 


EL 


AR 


EL 


12 

n 


(20) 

(21) 


The  unambiguous  range  interval  is 
R  =  c/(2*PRF) 

3 


(22) 


where  PRF  is  the  pulse  repetition  frequency.  The  number  of  range  bins  in 
an  unambiguous  range  interval  is  then 


where  Ar  is  the  range  bin  size  aid  is  determined  by  the  radar  video 
sampling  rate.  The  range  bin  index  for  range  Rmin  is 

„m.n  .  !NT  l<M0D[Rmjn,  R,  ])«/«,)  (24) 

where  M00(x,y)  is  equal  to  the  remainder  of  the  division  x/y,  and  INT(x) 
is  equal  to  the  largest  integer  whose  magnitude  does  not  exceed  x. 

The  minimum  range  contributing  to  the  m'th  range  bin 


R  ,  =  R  .  +  (m-m  .  j Ar  ;  m  >  m  . 

mO  min  ^  min  -  min 

(25) 

R  =  R  .  +  (M+m-m  .  W ;  m  <  m  . 

mO  min  inin  min 

(26) 

We  can  now  determine  if  the  range  bin 

of  interest 

( i.e.. 

the 

bin  index  rn)  is  illuminated  by  the 

main  bean. 

If  P-mO 

i  s 

qroater  than  Rn,JX,  the  main  beam  does  not  return  any  :lutter  in  bin 
number  n.  If  k^n  is  less  than  or  equal  to  there  will  be  one  or 

more  ambiguous  returns  in  bin  m.  If  the  number  of  ambiguous  range  returns 
is  small,  this  method  does  not  transition  smoothly  with  respect  to  bins 
just  inside  and  outside  the  beam  footprint.  If  desired,  the  transition 
may  be  made  more  accurate  by  considering  the  fraction  of  the  transmitted 
pulse  which  overlaps  the  footprint. 

The  number  of  contributing  ambiguous  range  returns  is 

K  =  I  NT { AR  /R  )  .  (27) 

'  FL  r  ' 


The  range,  grazing  angle,  and  elevation  angle  to  the  k'th  ambiguous  area 
that  contributes  to  the  m'th  range  bin  are  given  by 


The  illuminated  area  for  the  k'th  ambiguous  ranqe  and  m'th  ranqe  bin  is 
given  by 


\k  *  SLjAZjCSC!*^)  ;  tan  (*mk)  >  (31) 


c  T, 


^  3^ink 


AZ,  sec(4>  )  ;  tan  (♦.)<  -- (3?) 


c  t./? 


where  -  time  duration  of  the  received  compressed  pulse. 


The  mean  RCS  for  this  area  is 


RCS  =  o0(  +  ,  )  *  A 
rnk  mk  mk 


(33) 


where  o0  is  evaluated  usinq  Equation  1  or  2.  The  mean  power  for  the  m’th 
range  bin  is  then  qiven  by  the  son  of  the  power  from  the  contributing 
areas 


P  X2  r  K-l 


(4n) Jxc  k  =0 


G  (0,  Fl  -  FLo)  RC13  .  /R  ^ 
mk  ’  n<  mk 


(34) 


where 


□ 

t 


transmitted  power 
antenna  2 -way  power  qain. 


time  duration  of  the  transmitted  nulse. 


We  wish  to  remark  at  this  time  that  the  above  algorithms  are 
valid  in  a  statistical  sense  only.  While  we  have  decided  deterministical¬ 
ly  whether  a  particular  range  bin  is  illuminated  by  the  main  beam,  we  have 
neglected  to  account  for  certain  effects  such  as  atmospheric  refraction 
which  would  cause  an  offset  in  range  bin  position  relative  to  the  main 
beam  footprint.  In  support  of  this  approach,  we  note  that  one  range  bin 
has  no  greater  importance  than  any  other  as  far  as  radar  tracking  perfor¬ 
mance  is  concerned;  we  are  concerned  only  with  the  statistical  effect  on 
any  one  bin. 

4.2  MAIN  BEAM  DOPPLER  CENTER  FREQUENCY 

(neglecting  the  earth's  rotation) 

Without  implementation  of  QPSP,  the  mean  doppler  frequency  of 
the  main  beam  clutter  spectrum  for  the  m'th  range  bin  is  given  by 

2v 

f  =  — --  cos (EL  )  cos(AZ  )  (35) 

m  X  m '  v  o 

where 


El  =  -  ’  EL  . 

m  K  k=0  mk 


is  the  average  elevation  angle  of  the  ambiguous  returns  for  the  m'th  range 
bin.  Equation  36  gives  the  elevation  angle  as  the  average  over  the 

returns  from  the  K  ambiguous  ground  areas  that  arrive  simultaneously  to 

contribute  to  the  m'th  range  bin. 

With  QPSP  the  mean  frequency  of  the  ambiguous  returns  depends 

upon  the  pulse  to  which  the  inverse  phase  processing  is  matched.  For 

modeling  purposes  the  exact  value  is  not  important  and  as  a  convention  we 


fm  "  -  cos:ELm0;  cos;A70!  . 

X 


4.3  MAIN  BEAM  OOPPLcR  SPREAD 

(neglecting  the  earth’s  rotation) 

The  doppler  standard  deviation  due  to  internal  motion  is  given 
by  Equation  10  with  ov  nominally  4  m/sec. 

Doppler  standard  deviations  due  to  platform  motion  are  given  by 
.84  v 

oftZ  = - i  cos(Elo)  sin(A2o)  AZ,  (38) 


.84  v 

°EL  =  - ~  sin(EL0)  cos(AZo)  AEL  (39) 

X 

where 

AEL  =  Maximum  ( | ELrnK“ELml  |  ’  C^ctan ( ♦)/2R<. )  (40) 

without  QPSP,  and 

AEL  =  Minimum  (EL3,  Cretan  ( )  /  2R^  )  (41) 

with  QPSP. 

The  net  variance  of  the  spectrum  is  obtained  by  summing  the 
variances  due  to  platform  motion  and  to  internal  motion  (Equations  7  and 
11). 

4.4  DOPPLER  EFFECTS  INCLUDING  THE  EARTH'S  ROTATION 

Referring  to  Figures  6a  and  6b,  it  is  assumed  that  the  SBR 
location  is  given  in  geocentric  coordinates  by 
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Figure  6a.  SBR  and  earth  rotational  velocities  for  doppler  calculation 


Figure  6b.  SBR  LOS  In  the  tangent-plane  coordinate  systen. 


yV mm  ir»mW-rt  ■V-arVa-* 
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(-H  <  01  <  n) 


=  SBR  longitude 
4i  =  SBR  colatitude 
h  =  SBR  altitude 

and  that  the  azimuth,  elevation  angles  of  the  SBR  LOS  are  given  In  the  x, 
y,  z  tangent  plane  coordinate  system  by 

AZ  =  LOS  azimuth 

EL  =  LOS  elevation  (negative  when  looking  down)  . 

The  tangent -plane  coordinate  system  is  a  cartesian  coordinate 
system  with  x  directed  to  the  east,  y  to  the  north,  and  z  vertically. 

Vie  must  first  find  the  geocentric  coordinates  (a2.^2)  of  the  LOS 
intersection  with  the  earth's  surface. 

The  location  of  the  intersection  in  the  x',  y',  z'  cartesian 
coordinate  system  is  given  by 

x 2 '  3  {R  sin(ci)  sin(EL)  -  R  tan(Oi)  cos(EL)  sln(AZ) 

-  R  cosUi)  cos  (EL)  cos  (AZ)  +  (h+r0)  sin(c1)) 

/{cos(di)  <  sin(0! )  tan(0i)}  (42) 

y2‘  3  R  cos(EL)  sin(AZ)/cos(Oi)  +  tan(91)  x2‘  (43) 

for  J  3  v |  n/4  or  |3jj  >  3*/4  (case  a)  and  by 

y i'  =  lR  s i n L C i )  sin(EL)  +  R  cot(Ol)  cos(EL)  sin(AZ) 

-  R  cosUi)  cos(EL)  cos(AZ)  +  (h+r0)  slnicjl 

/  (sin  (  )  *■  cos  (  )  cot  (  Oj  j  [  (44) 

x2‘  =  -  R  cos  (EL)  sin(AZ)/sin(Oj )  +  cot('3L)  y2‘  (45) 
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for  n/4  <  j  © 1 1  <  3*/4  (case  b)  where  R  is  the  slant  range  found 
through  Equation  17  and 

22-  *  Iro2-  l<2')2-  (y2')2|1/2  .  (46) 

The  geocentric  coordinates  of  the  intersection  point  (92»42)  ace 

then 


02  *  tan- 1 (y 2‘ /x  2'  j 

(47) 

c.2  •  cos- l(.z2'/r0)  . 

(48) 

4.5  MAIN  8EAH  DOPPLER  CENTER  FREQUENCY 

The  doppler  frequency  of  a  point  on  the  earth's  surface  is 


fm  =  \  IVV*  (49) 

where 

vs  =  satellite  velocity  vector 
ve  =  earth's  surface  rotational  velocity 
R  =  unit  vector  along  the  LOS. 


The  satellite  velocity  vector  in  the  x,  y,  z  coordinate  system  is  assumed 
to  be  known  and  is  given  by 


v  =  v  „x  +  v  v  + 
s  sx  sy- 


vz 


R  is  found  to  be 


(50) 


R  =  sin(AZ)  cos(EL)x  +  cos(AZ)  cos(El)x  +-  sin(EL)z 


(51) 


(52) 


In  the  geocentric  coordinate  system 

_  A 

ve  “  “er°  s1n(52)0 

where  we  Is  the  earth's  radian  frequency  of  rotation. 


By  transforming  ve  to  the  x,  y,  z  coordinate  system  th_  dot  product  In 
Equation  49  may  be  evaluated  to  obtain 


f 

m 


>er0  sinU2)  cos ( 82- 9i ) )  sin(AZ)  cos(Et) 


+  ^vSy"“er°  sinUz)  cosUi)  sin(02-0i))  cos(AZ)  cos(EL) 

+  ^Vsz+“er°  s1n^2)  sinUi)  s1n(e2-0i))  sin  (EL ) }  . 

(53) 


4.6  MAIN  BEAM  DOPPER  SPREAO 


Doppler  spread  due  to  beamwldth  effects  may  be  found  by  taking 
the  deriative  of  Equation  49  with  respect  to  AZ  and  EL. 


°AZ 

=  .42 

3f  /3AZ 
m 

az3 

(54) 

°EL 

•  .42 

3f  /BEL 
m 

AEL 

(55) 

where  AEL  is 

gi  ven 

by  Equation 

40,  without  QPSP  or  by  41 

if  QPSP  is 

implemented. 

Referri ng 

to  Equations 

53  and  42  through  48, 

fm  depends 

implicitly  on 

AZ,  EL  through  92, 

t2  as  well  as  explicitly. 

Expressing 

fm  in  the  form 

f  =  -  (A  sin(AZ)  cos(EL)  +  B  cos(AZ)  cos(EL)  +  C  sin(EL)}  (56) 
m  X 
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where 


The  doppler  spread  due  to  azimuth  beamwidth  is 

°AZ  =  — 4— --  jA  cos(AZ)  ccs(EL)  -  B  sin(AZ)  cos(EL)  +  sin(AZ)  cos(El) 

TV0  cos^2)  cosio^e!)  1^2  +  w  sin(;  j  sin(e2-e1)  ^2] 

L  3AZ  e  3AZ  J 

-  cos(AZ)  cos( EL)  •  L)er 0  cos(c2)  cos(Cl)  sin(e2-0!)  lil 
L  3AZ 

+  7°  sin(;2)  costCi)  cos (V9!)  +  sin(EL) 

3AZ  J 

•  U  r0  cos  (c2  J  sin(cj)  s  in  (e2-  0 1 )  if* 

L  3AZ 

+  w  rQ  sin  (c2 )  sin(cj)  cos(o2-9i)  l!*]  (  • 

3AZJ  I 

The  doppler  spread  due  to  elevation  beamwidth  is 


(60) 


EL 


~--A--L-  J'A  sin(AZ)  sin(EL)  -  B  cos(AZ)  sin(EL)  + 

"  I 


C  cos(EL) 


+  sin(AZ)  cos(EL) 


V°  cos(c2)  cos;o2-0ij  i5? 

3EL 


4  Vo  s’nU2)  sin;  62-9,  j 

3FL  J 


-  cos(AZ)  cos(EL) 


f.  •••  ' 
\  -  .  * 
l.V.  '  ,V 


R- 


*.  »  »* 

*  • „  * 

•*  * •  •*.  %# ■*< 
p .  r. " 


A  =  v  - 

SX 

V  0 

sin (c2)  cos(e2-6!) 

(57) 

*W-' 

InTis!.^ 

B  =  v  - 
sy 

Vo 

sin(;2)  cos(  Ci)s  in  (e2-  ej 

(58) 

C  =  v  + 
sz 

V  0 

sinU2)  sin(Cj)  sin(e2-el) 

(59) 

v.v.v  % 
• .  • .  * 

t  .  •  . 

’  v.  .  •  J 


:*\*>-*j 


r  '  • 


r*  '*1 
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•‘vS . 


•  ,  «  , 
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»*• 
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•  <*>  r  o  COsU2)  COS  C«x  )  S1"n(e2~ei)  “ 

le  3EL 

Aft  1 

+  a)  r0  sin(c2)  cos  ( C x )  cos(e2-0i) _ 2  +  sin(EL) 

e  3ELJ 

*  U  r  q  cos  (c2)  5in(cJ  $10(62-6!)  iii 

L  e  3EL 


3EL 

+  ui  r0  sin(c2)  s in ( c ! )  cos(62-0!)  _^2 
e  3ELJ 


!  • 


(61) 


v'O 

-.i 


1 

1 


w. 

.*  V.  .* 
1^.-’ 


L*. 


»  ,  *  %  *  7 


The  partial  derivatives  of  02,  c2  ma>  *5e  found  from  Equations  47  and  48 

3*iL  y2'  li’) 

aAZ  anz  (. 

TTF) 


302  .  _ 1 

3AZ  1  +  (y2'/x2'  )2 


30, 

and  similarly  for  — 6 
3EL 


.aJL2'+  y  iii' 


x,'  — i  +  y  2‘ 

3C2  _  3AZ _ 3AZ^ 

3AZ 


^0*2'  A-(z2'  /r0)2 


(62) 


(61) 


t- 


•  *.  *  •  ■  1 


and  similarly  for 


3  4  2 
3EL 


The  remaining  chore  is  to  solve  for  the  partial  derivatives  of 
x 2*  *  y 2 ’  using  Equations  42  throuqh  46.  First  we  obtain  the  derivatives 
with  respect  to  azimuth.  For  case  a  we  have 


%  ’l* 


3*  2  ’ 
3AZ 


-  (-T  t  an  i.  0 ;  j  cns(EL)  cos(AZ)  +  R  cosLC.)  cos(El.)  sin(AZ)  } 


/  {cos  (ej  )  *■  sim'o.j  tan,  e  x  1 } 


(64) 
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=  R  cos(EL)  cos(AZ)/cos (et)  +  tanei  - 

3AZ  3AZ 


(65) 


and  for  case  b 


3y2' 


3AZ 


=  {R  cot (0 i )  cos(EL)  cos(AZ)  +  R  cos(ci)  cos(EL)  sin(AZ)} 


/(sin(ei)  ♦  cos(Qi)  001(0!)} 


(66) 


3x  1  3y  1 

— i-  s  -R  cos(EL)  cos(AZ)/sin  (0j )  +  cot(0i)  — —  . 

3AZ  3AZ 


(67) 


For  the  derivative  with  respect  to  elevation  we  first  need 


3R 


3EL 


=  - (r  0+h  )  cos(EL)  -  (r0+h)2  sin(El)  cos(EL) 


//(r0+h)2  s in 2( EL)  -  (2rgh+h2)  . 


(68) 


The  elevation  derivatives  are  then  given  for  case  a  by 


3x2' 

3~EL~ 


(cos ( 0 1 )  +  sin(0t)  tan(0i)}-1 

•  [s in (c  1 )  s in( £L)  -  tan(0j)  cos(EL)  sin(AZ) 

(  3EL 

-  005(4!)  cos(EL)  cos(AZ)  ]  *  R[sin(ci)  cos(EL) 

+  tan(ej)  sin(EL)  sin(AZ)  +  cos  (d)  sin(EL)  cos(AZ)]| 

(69) 


1*jL  -  sec (ei)»  (3R/3EL.  cos(FL)  sin(AZ)  -  R  sin(EL)  sin(AZ)} 


and  for  case  b  by 


l 

I 

E 

i 


=  {s in (e x )  +  cos(qJ  cot(ei))-1 
3EL 

•j*L-  Isin(cj)  sin(EL)  ♦  cot^}  cos(EL)  sin(AZ) 

( 3EL 

-  cos (c^)  cos(EL)  cos(AZ)  ]  +  R[sin(ci)  cos(EL) 

-  cottej  sin(EL)  sin(AZ)  +  005(5!)  sin(EL)  cos(AZ)]j  . 

(71) 

=  esc  ( ®  1 )  •  — — —  cos(EL)  sin(AZ)  ♦  R  sin(EL)  sin(AZ)!  . 

3EL  (  3EL  ' 

(72) 

In  sunmary,  these  algorithms  may  be  used  tc  find  the  doppler 
spreads  by  first  calculating  3R/3EL  and  then  the  partial  derivatives  of 
x 2 1  and  y2‘  with  respect  to  AZ  and  EL.  With  these  quantities  the  partial 
derivatives  of  e2,  C2  with  respect  to  AZ,  EL  may  be  computed  and  then  used 
to  solve  for  the  doppler  spread  through  Equations  64  and  55. 


Practically  speaking,  it  is  probably  computationally  less 
burdensome  to  solve  for  the  doppler  spreads  by  differencing  Equation  53 
over  the  beamwidths,  although  this  approach  does  require  solving  for  the 
four  earth  intersection  points  which  define  the  mainbeam  footprint.  Before 
choosing  either  method,  the  accuracy  and  computational  efficiency  of  each 
should  be  compared.  For  low  altitude  orbits  the  simpler  equations  derived 
without  inclusion  of  the  earth's  rotational  effects  are  approximately 
correct  and  may  be  utilized. 


4.7 


QUADRATIC  PHASE  SHIFT  PROCESSING  (QPSP) 


The  purpose  of  QPSP  is  to  reduce  the  dopp'ier  spread  of  the  main 
beam  clutter  spectrun  by  compensating  for  the  shift  in  doppler  center 
frequencies  of  the  ambiguous  range  returns.  Modeling  the  doppler  spread 
with  or  without  QPSP  was  discussed  earlier  in  this  section.  We  now 
describe  the  basis  for  the  signal  processing  involved  in  QPSP. 


In  this  technique  the  phase  is  shifted  quadrat ical ly  in  one 
direction  as  a  function  of  pulse  number  on  transmit,  and  in  the  opposite 
direction  on  receive.  We  assume  that  a  phase  ^nift  modulation 


=  ep 


(73) 


is  implemented  on  a  pulse-to-pulse  basis  within  a  sequence  of  transmitted 
pulses  where  p  is  pulse  number.  On  receive,  the  range-gated  signal  may 
contain  main  beam  clutter  power  returned  from  one  or  more  of  the  trans¬ 
mitted  pulses.  When  there  is  more  than  one  contributing  pulse,  these 
multiple  and  simultaneously  received  pulses  are  called  ambiguous  range 
returns.  We  use  the  index  k  once  again  to  designate  the  ambiguous  range 
number,  with  k  increasing  with  increasing  range,  and  choose  a  time 
reference  such  that  the  sampling  time  when  the  p'th  pulse  is  the  zeroth 
ambiquous  return  is  pT.  As  an  example,  in  Figure  7  we  show  the  relation 
of  pulse  number  to  ambiguous  range  index  and  sampling  time  for  the  case 
that  there  are  five  contributing-  ambiguous  returns.  Then  the  phase  on  the 
k'th  ambiquous  return  at  time  pT  is  given  by 


$tk(pT,k)  =  c(P*k)2  •  (74) 

We  assume  that  a  quadratic  phase  shift  of  the  opposite  siqn  is  applied  on 
receive 

<t>r  { pT ,  k )  -  -t(p- Apl  2  (75) 
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Pulse  mmtoer  versus  sanpling  tine  and 
ambiguous  range  index. 


when?  the  offset  Ap  reflects  the  fact  that  the  inverse  phase  shift 
processing  in  general  may  not  be  matched  to  that  of  the  transmitted  phase 
modulation . 


is  then 


The  net  phase  modulation  on  the  k'th  ambiguous  return  at  time  pi 


‘t'net  ^pT *k  >  *  c(k?-Ap+2p  Ap-2pk  ) 


and  the  doppler  shift  of  the  k'th  ambiguous  return  is  then 

i  dVt(p!''k) 

Af  (k  )  = - - =  cAp/irT  -  ek/nT  .  (77) 

2n  dpi 


QPSP  is  seen  to  impose  a  frequency  shift  that  varies  linearly  with  ambigu¬ 
ous  range  index. 

The  doppler  shift  of  the  k'th  ambiguous  return  due  to  the  eleva¬ 
tion  angle  difference  is  given  by 


Af  (k)  =  sin  (EL0 )  cos(AZo)  AEL(k) 
CL  X 


where 


AEL(k)  -  R  k  tan(i(/0)/R 

3  C 


For  ease  of  interpretation  we  have  used  the  expression  for  Afc|_  neglecting 
the  earth's  rotational  effects.  For  the  frequency  shift  from  processing 
to  cancel  that  due  to  ambiguous  range  we  require 


Af  pi  (k )  =  -sk/fT 
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since  the  first  term  in  Equation  77  is  a  frequency  offset  which  Is 
constant  with  respect  to  ambiguous  range.  Using  Equations  78  and  79  to 
solve  Equation  77  for  e  we  obtain 


e  =  -2nv  R  T  sin(EL0)  cos(AZ0)  tan(^0)/xRr 

S  d  t 


If  we  had  used  the  expression  for  Af£|_  that  accounted  for  the 
earth's  rotational  effects  in  the  above  derivation,  e  would  still  have 
been  calculable  in  the  same  manner. 

In  summary,  QPSP  has  been  shown  to  be  an  effective  technique  for 
mitigating  the  doppler  shift  of  ambiguous  returns.  This  has  the  desirable 
result  of  increasing  the  clutter  free  doppler  bandwidth  available  to  a 
pulse  doppler  radar.  When  applied  in  an  optimum  sense  it  requires  imple¬ 
mentation  as  a  function  of  look  angle. 


SECTION  5 


SIOELOBE  CLUTTER  MODEL 


This  section  is  devoted  to  the  calculation  of  the  clutter  power 
spectrum  that  is  caused  by  returns  from  outside  the  main  beam  footprint. 
Because  of  an  SBR's  large  field  of  view,  the  clutter  signal  received 
through  the  antenna  sidelobes  will  contain  multiple  ambiguous  range 
returns.  The  illuminated  area  for  the  sidelobe  return  in  a  given  range 
bin  is  a  series  of  annular  rings  corresponding  to  the  pulse  length  limited 
areas  of  the  ambiguous  returns.  Mean  power  of  the  sidelobe  return  may  be 
taken  to  be  the  same  for  all  range  bins,  since  the  ambiguous  returns  for 
different  bins  will  sample  essentially  the  same  areal  portion  of  the 
earth's  surface. 

To  calculate  the  mean  power  we  must  sin  the  power  contributions 
from  the  ambiquous  returns.  It  is  desirable  to  estimate  the  range  at 
which  the  power  returned  becomes  negligible  so  that  the  number  of 
ambiguous  returns  to  be  summed  may  be  reduced  accordingly.  The  power  per 
ambiguous  return  falls  off  at  least  as  fast  as  the  reflectivity  function 
00.  By  referrinq  to  Figure  3  we  see  that  for  both  the  HW  and  LW  models 
o0  has  decreased  by  more  than  20  dB  f^om  its  peak  value  at  a  grazing  angle 
of  50  degrees.  The  range  at  which  this  occurs  is  given  by 

°50  =  -r0sin(50°)  +  l/^r0^in2(50’)  +  (2r0h+hz)  .  (82) 

Thus,  we  need  only  consider  returns  for  ranges  less  than  R  .  As  an  Ex¬ 
ample,  for  a  low  altitude  (1400nm)  orbit  and  a  PRF  less  than  10  KHz,  fewer 
than  17  ambiguous  returns  are  important  for  calculation  of  the  required 


Figure  8.  Geometry  for  sidelobe  doppler  analysis 


10000.  20000.  30000. 

ORBIT  ALTITUDE  (KILOMETERS) 


*0000. 


Figure  9.  Sldelobe  clutter  doppler  shifts. 


grazing  angle  of  80  degrees  as  a  function  of  orbit  altitude.  Also  shown  on 
the  fiqure  is  the  one  sided  spectral  width  a  .  .  We  note  that  the  two 

sided  spectral  width  is  greater  than  10  KHz  for  orbit  altitudes  less  then 
-5,000  kilometers  (-2,700  nm) .  Since  we  are  primarily  interested  in  this 
range  of  orbits  and  for  PRFs  equal  to  about  10  KHz  and  smaller,  the 
sidelobe  return  may  be  modeled  as  a  flat  spectrum,  similar  to  additive 
white  noise. 

5.1  SIDELOBE  CLUTTER  POWER  IN  THE  M'TH  RANGE  BIN 


The  number  of  ambiguous  range  returns  considered  is 


‘si  ■  INT!<RS0-  h,/Ra| 


(86) 


where  R50  is  obtained  through  Equation  82. 


The  range,  qrazinq  angle,  and  illuminated  area  of  the  k'th 
ambiguous  return  are  given  by 


Rk  =  h  +  (k-1)  •  Ra 


4;  =  S  in 


.  I 


l!u  ♦  l.  j  -  "a.  ! 


—  ■  ■ - ) - i 

Rkl  2ro  2ro  I 


A,  =  n  c  x  R, 
k  c  k 


(87) 

(88) 

(89) 


The  mean  power  of  the  k'th  return  is 


V2GSLC’t  °0(V 


64  iV 


(90) 


where  Gt-L  is  the  average  two-way  antenna  sidelobe  gain. 


The  mean  power  in  any  range  bin  is  then  given  by 
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APPENDIX 


A  STATISTICAL  SIGNAL  GENERATION  TECHNIQUE 
FOR  REALIZING  THE  SA*»LE  CLUTTER  VOLTAGE 

(0.  L.  Knepp) 

In  this  appendix,  a  numerical  technique  is  described  to  generate 
realizations  or  sample  functions  of  the  video  clutter  voltage  received  by  a 
pulse  doppler  SBR.  The  goal  of  this  technique  is  to  generate  the  two  quad¬ 
rature  components  of  the  clutter  siqnal  as  received  at  a  single  range  gate 
for  all  the  pulses  of  a  coherent  pulse  train  or  coherent  dwell.  It  will  be 
shown  that  this  technique  requires  knowledge  of  the  mean  clutter  spectrum 
and  depends  upon  the  observation  that  the  clutter  voltage  has  a  Gaussian 
probability  distribution  with  uncorrelated  quadrature  components.  An  iden¬ 
tical  statistical  signal  technique  is  used  by  Knepp1  to  generate  realiza¬ 
tions  of  the  phase  for  a  multiple  phase  screen  propagation  simulation. 

Assume  that  there  are  N  pulses  oer  coherent  pulse  train  separ¬ 
ated  in  time  by  the  quantity  T  so  that  the  pulse  repetition  frequency  (PRF) 
is  1/T.  As  will  be  shown,  the  clutter  voltage  c(nT),  is  generated  from 
initial  knowledge  of  its  mean  power  spectrum.  In  continuous  notation,  the 
clutter  voltage  may  be  written  as  the  Fourier  transform 

c(t)  =  jf  C(w)e1u,td<D  ( A-l) 


1.  "Multiple  Phase-Screen  Calculation  of  the  Temporal  Behavior  of  Stochas¬ 
tic  Waves,"  Proc.  IEEE,  Vol.  ?1,  No.  6,  pp.  722-737,  June  1983. 
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In  the  discrete*  case.  Equation  A-l  is  written  as 


c(nT)  =  *1  C(fflHel2w/N  Au  n=0 . N-l  (A-2) 

ma0 

Now  if  the  Fourier  transform  of  the  clutter  spectrun  was  available,  a 
clutter  voltage  realization  could  be  easily  generated  by  using  Equation 
A-2.  For  the  moment  let  us  choose  as  the  Fourier  transform  the  quantity 

C(mAu>)  =  rm[C(mAw)NT/2it]1/  2  (A-3) 

where  C(mAu>)  represents  discrete  values  of  the  known  desired  mean  clutter 
power  spectrum.  In  the  following  it  is  proven  that  this  choice  is 
correct.  Here  Au>  =  2*/NT,  NT  is  the  time  duration  of  the  dwell,  and 
therefore  TAm  =  2"/N.  r^  is  a  complex  nunber  given  as  the  sum  of  two  inde¬ 
pendent  Gaussian  random  variables  with  zero  mean  and  variances  of  unity. 


Successive  values  of  glfT)  and  may  be  obtained  numerically  by  sampling 
from  a  pseudo-random  sequence  of  numbers  with  a  Gaussian  distribution.  The 
factor  of  Kii  is  included  so  that 


<'rmrn>  =  ®mn  ( A-  5) 

wh-r°  'Wi  ’s  the  Kroneker  delta  function.  It  is  apparent  that  with  the 
above  choice  for  rm,  c(nT)  is  the  sun  of  a  sequence  of  Gaussian  variates 


*  In  this  appendix  all  discrete  sums  are  taken  over  a  range  of  the  index 
from  0  to  N-l. 
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and  thus  its  real  and  imaginary  parts  (i.e.,  its  quadrature  components) 
both  have  a  Gaussian  or  normal  probability  distribution. 


So  that  Equation  A-8  becomes 


Bc(kT)  =  l  C(mAu))ei2imlk/N  IrJ 2  to  (A-10) 

m  1  ' 

A  comparison  of  the  above  equation  to  the  continuous  relationship  between 
power  spectrum  and  autocorrelation  function  as  given  by 

Bc(t)  =_£  C( w)eia,tdu>  (A- 11) 

shows  that  the  power  spectrum  of  the  numerically  generated  clutter  voltage 

is  C(mAai)  |rm|  2.  It  is  apparent  that  different  values  of  the  index  m  corre¬ 
spond  to  different  radian  frequency  components  of  the  power  spectrum. 

Since  |rm|2  is  the  sum  of  the  squares  of  two  Gaussian  variates,  each  of  the 
Fourier  components  of  the  power  spectrum  of  an  individual  clutter  realiza¬ 
tion  is  a  chi-squared  variate  with  two  degrees  of  freedom  and  a  mean  value 
of  C(mAu>).  Thus  for  any  given  clutter  realization,  the  power  spectrum  will 
not,  in  qeneral,  be  identical  to  the  desired  spectrum.  However,  the  aver¬ 
age  power  spectrun  of  many  such  clutter  voltage  real  izaticwis  may  be  ob¬ 
tained  by  taking  the  expected  value  of  Equation  A-10.  Since  <|rm|2>  =  1 

AVG[Bc(kT)  ]  =  ^  C(mAu>)el2im,k/NAu>  (A-12) 

which  may  be  compared  to  Equation  A-ll.  Therefore  the  average  spectrum 
obtained  from  many  clutter  voltage  realizations  is  the  desired  power 
spectrum,  C(mAu)),  and  Equation  A-3  if  correct. 
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Joint  C'uise  Missiles  Project  Ofc 

ATTN: 

JCMG-707 

Naval  Air  Systems  Ccr.Tiiand 

ATTN: 

PMA  271 

Naval  Electronic  Systems  Command 

ATTN: 

PME  106-4,  S.  Kearney 

ATTN: 

Code  S01A 

ATTN: 

PME  117-211 .  C.  nruger 

ATTN : 

Code  3101,  T.  Hughes 

ATU. 

PME  - 106 ,  F.  ijiederich 

ATTN: 

PME- 1 17-2013,  G.  Burnhart 

ATTN: 

PME  117-20 

Naval  Intelligence  Support  Ctr 

ATTN: 

K I  SC  *30 

Naval  Ocean  Systems  Center 

ATTN: 

Code  532 

ATTN: 

Code  S322,  M.  Paulson 

ATTN: 

Code  5323,  J.  Ferguson 

Naval  Research  Laboratory 

ATTN: 

Code  -720,  J.  Davis 

ATTN: 

Code  4700,  s.  jssaicv. 

ATTN: 

Code  4 1  Go ,  L.  : zjstewic l 

ATTN: 

Code  /SOD,  b.  /.aid 

AT  Til : 

Code  6700 

r-sT  7N : 

code  4/00 

AT  :N  ; 

oo  ce  ‘  7  JO 

AH  ?l : 

Code  7HSO,  J.  aj-'Cnii. 

ATT.’i: 

Cooe  4 1  <j  7 

Naval  S^jce  Surveil  lanct?  System  -■ 

ATTN : 

J.  Burton 

Naval  Surface 

,.0.1  cons  Center 

ATTN: 

Code  FJ1 

Naval  Tr.-lec-j.-i 

l Cd  t  ions  r  .  i-  J ■  ■  cl 

ATTN: 

Cooe  341 

C'fc  of  tnc  uc 

pc’  -/  Cnic*  of  No/al  Ops 

AT  IN 

No  Nh  1  j 

mi  i  : i : 

NO?  ?J1N 

ATlfJ: 

.,'jP  ;-3v  ,  >  t '  d  t  c  ■  j  1  :•  mNi) 

DEPARTMENT  OF  THE  NAVY  (Conti nued) 

Office  of  Naval  Research 

ATTN:  Cod.-  414,  G.  Joiner 
ATTN:  Code  412,  W.  Condel 1 

Strategic  Systems  Project  Office 

ATTN:  NSP-4J,  Tech  Library 
ATTN:  NSP-21 41 
ATTN:  NSP-2722 

Theater  Nuclear  warfare  Proj  Office 
ATTN:  PM-23,  0.  Smith 

DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Geophysics  laboratory 
ATTN:  LYD,  K.  Champion 
ATTN:  OPR- 1 
ATTN:  R.  Sabcoci. 

ATTN:  CA.  A.  Stair 
ATTN:  OPR,  H.  Gardiner 
ATTN:  PHY.  J.  3uchau 
ATTN:  R.  O'Neil 

Air  Force  Satellite  Ctrl  Facility 

Uet  3,  AWS 

ATTN:  WE 

Air  Force  Space  Technology  Ctr 
ATTN-  YH 

Air  Force  Technical  Applications  Ctr 
A1TN:  TN 

Air  Force  -.teapons  Laboratory 
ATTN:  SUL 
ATTN:  NTS 
ATTN:  NTCA 

Air  Force  Wright  Aeronautical  Lat/AAAO 
ATTN:  VI.  Hunt 
ATTN:  A.  Johnson 

Air  Log i  sties  Corr.and 
ATTN:  OO-AlC/HK 

Ai-  University  Library 
ATTN:  Al-L -L  SC 

Assistant  thief  of  Staff 

Studies  &  Ana  1 ys ’s 

ATTN:  AI/SASC,  C.  Rignto'eyer 

ballistic  Missile  Office/DAA 
ATTN:  Sit,  U.  tv, an 
ATTN:  ENSN 

ATTN:  EilS'i,  vi .  Wilson 

Deputy  Cnief  of  Staff 

Research.  Dev  el  opr  r-nt .  i,  Acg 
ATTN:  AFPD-VI 

ATTN:  ArSJS  ,  Space  Sys  i  C.3  Ji 

Deputy  Chief  of  staff 

Plans  dii-J  Gi.-ei  iti/ns 
ATTN:  AF'OtCD 
ATTN:  Ar>:-.T 
ATTN-  A  r 1 0 ' . 
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nEPARTHENT  OF  THE  AIR  FORCE  (Continued) 

Electronic  Systems  Div 

ATTN:  ESO/SCTE,  J-  Clark 

Electronic  Systems  Division 

ATTN:  SCS-2,  G.  Vinkels 
ATTN:  SCS-1E 

foreign  Technology  Division 

ATTN:  TQTD,  8.  Ballard 
ATTN:  Nits,  Library 

Rome  Air  Development  Center 
ATTN:  OCS,  V.  Coyne 
ATTN:  OCSA ,  R.  Schneible 
ATTN:  TSLO 

Rome  Air  Development  Center 

ATTN:  EEP,  J.  Rasmussen 
ATTN  EE PS,  P.  Kossey 

Space  Command 

ATTN:  DC,  T.  Long 
ATTN:  XPSD,  J.  Hatlel id 

Strategic  Air  Command 
ATTN:  AOWA 
ATTN:  OCX 
ATTN:  XPQ 
ATTN:  XPFC 
ATTN:  OCZ 

ATTN:  MR  I/ST  INFO  Library 
ATTN:  XPFS 

DEPARTMENT  OF  ENERGT 

Department  of  Energy 
GTN 

ATTN:  OP-233 

OTHER  GOVERNMENT  AGENCIES. 

Central  Intel  1 igence  Agency 
ATTN:  OSWR/NED 

ATTN:  OSWR/SSO  for  K.  Feuerpfetl 

Department  of  Commerce 
National  Bureau  of  Standards 

ATTN:  Sec  Ofc  for  R.  Moore 


DEPARTMENT  OF  ENERGY  CONTRACTORS 

EGSG,  Inc 

ATTN:  J.  Colvin 
ATTN:  0.  Wright 

University  of  California 

Lawrence  Livermore  National  Lab 
ATTN:  L-31 ,  R.  Hager 
ATTN.  Tech  Info  Dept  Library 

Los  Alamos  National  Laboratory 
ATTN:  T.  Kunkle,  ESS -5 
ATTN:  0.  Simons 
ATTN:  MS  664,  J.  Zinn 
ATTN:  MS  670,  J.  Hopkins 
ATTN:  R.  Jeffries 
ATTN:  J.  Wolcott 
ATTN:  P.  Keaton 

Sandia  Nation.  I  Laboratories 
ATTN:  T.  Cook 
ATTi  .  B.  Murphey 

Sandia  National  Laboratories 

ATTN:  Org  1250,  U.  Brown 
ATTN:  Org  4231,  T.  Wright 
ATTN:  Space  Project  01  v 
ATTN:  0.  Thornbrouoh 
ATTN:  Tech  Library  3141 
ATTN:  C.  Oahlgren 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corp 

ATTN:  1.  Garfunkel 
ATTN:  0.  Olsen 
ATTN:  J.  Kluck 
ATTN;  0.  Whelan 
ATTN:  J.  Straus 
ATTN:  R.  Slaughter 
ATTN:  K.  Cho 
ATTN:  V.  Josephson 
ATTN:  T.  Salmi 

Aerospace  Corp 

ATTN:  S.  HcWaters 

Analytical  Systems  Engineering  Corp 
ATTH:  Radio  Sciences 
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Department  of  Comnerce 

National  Oceanic  6  Atmospheric  Admin 


Institute  for  Tel ecoiifnun icat ions  Sciences 
National  Telecommunications  &  Info  Admin 
ATTN :  A .  Jean 
ATTN.  L.  Berry 
ATTN:  W.  Utlaut 

Department  of  State 

Office  of  International  Security  Policy 
bureau  of  Politico  Military  Affjirs 
ATTN:  PM/STM 

NATO 

NATO  School,  SHAPE 

ATTN:  US  Documents  Officer 


DEPARTMENT  Of  DEFENSE  CONTRACTORS  (Continued 

University  of  California  at  San  Diego 
ATTN:  H.  8ooker 

Charles  Stark  Draper  Lab.  Inc 
ATTN:  A.  Tetewskl 
ATTN:  0.  Cox 
ATTN:  J.  Gilmore 

Communications  Satellite  Corp 
ATTN:  0.  Fang 
ATTN:  G.  Hyde 

Computer  Sciences  Corp 

ATTN:  T.  ElsenSarth 

Cornell  University 

ATTN:  0.  Farley,  Jr 
ATTN:  H.  Kelly 

E-Systems ,  Inc 

ATTN:  R.  Berezdivin 

Electrospace  Systems,  Inc 
ATTN:  P,  Phillips 
ATTN:  H.  Logs  ton 

EOS  Technologies,  Inc 

ATTN:  W.  Lelevier 
ATTN:  B.  Gabbard 

General  Electric  Co 

ATTN:  R.  Junpr 
ATTN:  A.  Steinmayer 
ATTN:  C.  Zierdt 

General  Electric  Co 

ATTN:  G.  Hillman 

General  Research  Corp 

ATTN:  B.  Bennett 
ATT?  R.  Rein 
AHN:  R  Williams 

GEO  Centers,  Inc 

ATTN:  E.  Ha rram 

GTE  Comnunl cat  Ions  Products  Corp 
ATTN:  R.  Stelnhoff 

GTE  Conmunlcatlons  Products  Corp 
ATTN:  J.  Concordia 
ATTN:  1.  Kohl  berg 

Harris  Corp 

ATTN:  E.  Knlck 

Honeywel 1 ,  Inc 

ATTN:  G.  Terry,  Avionics  Cept 
AT  1 N :  A.  Kearns  MS924-3 

Horizons  Technology,  Inc 
ATTN:  R,  Kruger 

HSS,  Inc 

ATTN:  0.  Hansen 
IBM  Corp 

ATTN:  H.  Ulander 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued 

Institute  for  Defense  Analyses 
ATTN:  H.  Gates 
ATTN:  J.  Aein 
ATTN:  E.  Bauer 
AHN:  H.  Wolf  hard 

ITT  Corp 

ATTN:  Tech  Library 
ITT  Corp 

ATTN:  G.  Wetmore 

JATCQR 

ATTN:  J.  Sperling 

JAVCOR 

AHN:  H  Oickinson 

Johns  Hopkins' University 
ATTN:  C.  Meng 
ATTN:  K.  Potockl 
ATTN:  J.  Phillips 
ATTN:  T.  Evans 
ATTN:  J.  Newland 
ATTN:  P.  Komlske 

Kaman  Sciences  Corp 
Systems  Directorate 

ATTN:  E.  Conrad 

Kaman  Tempo 

ATTN:  B,  Gamblll 
ATTN:  CAS  I AC 
ATTN:  w.  Schuleter 
ATTN:  W.  McNamara 

Kaman  Tempo 

ATTN:  DAS [AC 

Litton  Systems,  Inc 

ATTN:  B.  Zlrnier 

Lockheed  Missiles  4  Space  Co,  Inc 
ATTN:  R.  Sears 
ATTN:  J.  Kumer 

Lockheed  Missiles  4  Space  Co,  Inc 
ATTN:  Dept  60-12 
2  cy  AT1N  D.  Churchill 

M.  1 .  T  ,  Lincoln  Lab 

ATTN:  N.  Doherty 
ATTN:  V.  Vltto 
ATTN:  0.  Towle 

M/A  Com  Ltnkabit  Inc 

ATTN:  A,  Viterbi 
ATTN :  H .  Van  T  rees 
ATTN:  I.  Jacobs 

Magnavox  Govt  4  Indus  Electronics  Co 
ATTN:  G.  White 

Maxim  Technologies,  Inc 
ATTN:  E.  Tsui 
ATTN:  J,  Marshall 
ATTN:  R.  Horganstern 
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U€ PART HEfft  OF  DEFENSE  CONTRACTORS  (continued) 


DEPARTMENT  OF  PEFENSE  CONTRACTORS  (Continued) 


McDonnell  Douglas  Corp 

RAD  Associates 

ATTN:  Tech  Library  Services 

ATTN:  G.  Ganong 

ATTN:  R.  Halprin 

ATTN:  W.  Olson 

RAND  Corp 

Meteor  Communications  Corp 

ATTN:  E.  Bedrozlan 
ATTN:  P.  Davis 

ATTN:  R.  Leader 

ATTN:  C.  Crain 

Mission  Research  Corp 

P.and  Corp 

ATTN:  R.  Blgoni 

ATTN:  B.  Bennett 

ATTN:  G.  McCartor 

ATTN:  F.  Gulgllano 

Riverside  Research  Institute 

ATTN:  F.  Fajen 

ATTN:  v.  Trapani 

ATTN:  D.  Knepp 

ATTN:  Tech  Library 

Rockwell  International  Corp 

ATTN:  R.  Kllb 

ATTN:  R.  Buckner 

ATTN:  R.  Bogusch 

ATTN:  R.  Hendrick 

Rockwell  International  Corp 

ATTN:  S.  GutSChe 

ATTN:  S.  Qulllcl 

ATTN:  R.  Dana 

2  cy  ATTN:  C.  Lauer 

Science  Appl (cations,  Inc 

5  cy  ATTN:  Document  Control 

ATTN:  D.  Sachs 

Mitre  Coro 

ATTN:  C.  Smith 

ATTN:  D.  Hamlin 

ATTN:  G.  Harding 

ATTN:  L.  LinsOn 

ATTN:  A.  Kymnel 

ATTN:  E.  Straker 

ATTN:  MS  J104,  M.  Oresp 

ATTN:  C.  Callahan 

Science  Applications,  Inc 

Mitre  Corp 

ATTN:  W.  Hall 

ATTN:  J.  Cockayne 

Science  Applications,  Inc 

ATTN:  J.  Wheeler 

ATTN:  M.  Cross 

ATTN:  V.  Foster 

ATTN:  M.  Horrocks 

SRI  International 

Pacific-Sierra  Research  Corp 

ATTN:  A.  Burns 

ATTN:  G.  Price 

ATTN:  F.  Thomas 

ATTN:  R.  Tsunoda 

ATTN:  H.  Brode,  Chairman  SAGE 

ATTN:  J.  Vickrey 

ATTN.  E.  Field,  Jr 

ATTN:  V.  Gonzales 

Pennsylvania  State  University 

ATTN:  D.  He  11  son 
ATTN:  J.  Petrtckes 

ATTN:  Ionospheric  Researcn  Lab 

ATTN:  M.  Baron 

Photometries,  Inc 

ATTN:  R.  Livingston 
ATTN'  D  McDaniels 

ATTN:  I.  Kofsky 

ATTN:  W.  Chesnut 

Physical  Dynamics,  Inc 

ATTN:  G.  Smith 

ATTN:  C.  Rlno 

ATTN:  J.  Secan 

ATTN:  W.  Jaye 

ATTN:  E.  Fremcuw 

ATTN:  R.  Leadabrand 

Physical  Research,  Inc 

Stewart  Radtance  Laboratory 

ATTN:  R.  Dellberls 

ATTN:  J.  ulwlck 

ATTN:  J.  Devore 

ATTN:  T  Stephens 

Swerllng,  Manasse  S  Smith,  1 

ATTN:  J.  Thompson 

ATTN:  R.  Hanasse 

RtO  Associates 
ATTN 
ATTN 
ATTN 
ATTN 
ATTN 
ATTN 
ATTN 
ATTN 
ATTN 


Turco 

Gantsweg 

Grelf Inger 

Gilmore 

Ory 

Steyr 

Karzas 

Wrignt 

l,a.i  s 


Technology  International  Corp 
ATTN:  W.  Boqulst 

Toyon  Research  Corp 
ATTN:  J.  Ise 
ATTN.  J.  Garbarlno 

Utah  State  University 


R&D  Associates 

ATTN:  B.  Yoon 


ATTN 

ATTN 

ATTN 

ATTN 


A  Steed 
D  Burt 

K.  Baker.  Uir  Atmos  t>  Space  Sci 

L.  Jensen,  Elec  Eng  Dept 
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OEPARTMENT  OF  OEfENSE  CONTRACTORS  (Continued) 

Visidyne,  Inc 

ATTN:  j.  Carpenter 
ATTN:  W.  Reidy 
ATTN:  0.  Shepard 
ATTN.  C.  Humphrey 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

TRW  Electronics  l  Defense  Sector 
ATTN:  R.  Plebuch 

Boe(nq  Co 

ATTN;  MS  8K-8S,  Dr  S.  Tashird 
ATTN:  MS/87-63,  0.  Clauson 
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